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ABSTRACT 

We consider accretion onto newborn black holes (BHs) following the collapse of rotating massive 
stellar cores, at the threshold where a centrifugally supported disk gives way to nearly radial inflow for 
" low angular momentum. For realistic initial conditions taken from pre-supernova (pre-SN) evolution 

. calculations, the densities and temperatures involved require the use of a detailed equation of state 

' and neutrino cooling processes, as well as a qualitative consideration of the effects of general relativity. 

^Sj I Through two-dimensional dynamical calculations we show how the energy release is affected by the 

. rotation rate and the strength of angular momentum transport, giving rise to qualitatively different 

' solutions in limits of high and low angular momentum, each being capable of powering a gamma-ray 

^) I burst (GRB). We explore the likelihood of producing Fe-group elements in the two regimes and suggest 

■ that while large and massive centrifugally supported disks are capable of driving strong outflows with 

' a possible SN-like signature, quasi-radial flows lack such a feature and may produce a GRB without 

I such an accompanying feature, as seen in GRB060505. 

^ ' Subject headings: accretion, accretion disks — hydrodynamics — gamma rays: bursts — supernovae 

O ■ 1- INTRODUCTION 

-^-^ . Gamma-ray bursts (G RBs) have remained as an outstanding problem in astrophysics since their discovery in 1967 
^ ' (|Klebesadel et al1ll973f ). Major advances in their understanding came through observational breakthroughs following 
BATSE, which provided an extensive all-sky catalog (Fish man & Mcegan 1995) and the identification of counter- 

J iarts at lower energ ies at high redshift (|Metzger et al.|[l997l ) through Beppo Sax ()van Paradiis et aI]|2000D and Swift 
Gehrels et al.l [20071 ). Isotropic equivalent energies range from 10*^ erg s~^ to 10^^ erg s~^, coming from sources at 
»^ ' redshifts as high as z ~ 6 that are randomly distributed over the sky ([Meegan et al.l [l"992l ). GRBs are commonly 
divided into short (SGRBs, less than 2 seconds in duration, typically at redshift 2; ~ 1 and frequently in host galaxies 
with low star formation rates (SFRs)) and long (LGRBs, longer than 2 seconds, residing at z ~ 1 — 5 in star forming 
galaxies) events. An additional dist inction is that those of the short variety tend to have harder spectra, hence the 
short-hard/long-soft denomination (jPezalav et al.lfr996( l. Recently there appears to be a substantial and, in some 
00 ' respects, puzzling overlap between the two populations, and it has been suggested that a third group or even a new 

classification is necessary (iGehrels et aLllMeHGaTYam et al.l[200l iDella Valle et al.l[2006bl: iFvnbo et"al][200l . 
QP ■ The energy release, short duration, and variability time scales associated with GRBs favor accretion onto com- 
pact objects, or magnetically powered events from such objects, as their ultimate source. Many models have been 
suggested over the years, a number of which were discarded once the global energy budget was fixed through the 
" ■ resolution of the distance scale to the sources. Amo ng the early proposa ls still considered are comp act binary m ergers 
K>N ' jLattimer & Schramm" 1974', '1976'; 'Paczvhski![i98|, Il991t lEichler et anil989t iNaravan ei al."l992^. magnetars (Us^ 
1992; Thompson 1994; Mcszaros & Recs 1997,) and the co llapse of massive stellar cores (|Wooslcy 1993), sec (Mcszarog 
2002t |Pira3l2004l: lNakaill2007l : iLee fc Ramirez-Ruizl[200l for reviews. 



o 



o 



In general, accretion onto a central object can proceed and release gravitational binding energy efficiently if the 
gas c an cool or otherwise lose its in ternal energy e.g., through advection, and thus move down in the potential 
well (ISalpeteil 11964 [ZerDovichlll964f ). At the usual astrophysical scales, this occurs through photons, and produces 
phenomena ranging from active galactic nuclei (AGN) to low mass X-ray binaries (LMXBs). If the accretion rate is 
too high though, the fluid may become so optically thick to the photons that they become trapped and are unable 
to cool the gas. At even higher accretion rates, a different mechanism enters the picture, namely, neutrino cooling. 
This regime is termed "hypercritical" , since it is far above the usual Eddington limit for photons, an d is relevant, for 
example, in post-supernova cores and the associated fall back, as was probably the case in SN1987A (|Chevalieilll989l : 
iHouck fcC hevalier 1991). It is this phenomenon that is now thought possible of generically powering GRBs. The 
question at the forefront of attention has been how to achieve these conditions in an astrophysical setting. 

There has been mounting evidence in the past few years l inking LGRBs at low redshift with ty p e Ic SNe, with spatia l 
and temporal coincidences for a number of events (Wooslc v fc Blooii3l2006t iKaneko et al.ll2008l ). iGalama et al.l (|1998l ) 
hnked for the first time a LGRB with a SN (GRB980425 with SN1998bw). After 27 days the optical spectrum of 
GRB980425 resembled that of a H-deficient type Ic SN. To date there have been at least 6 other LGRB/SN connections: 



XRF020903 llSod erberg eral][2005[) : GRB021211 with SN20021t (jPella Valle et all [20031): GRB030329 with SN2003dh 
(whose optical spe ctrum afterglow was strikingly simi lar to that of SN1998bw (jStanek et a l. 2003)); GRB031 203 wit h 
SN20031W (Malcsani ct al."2004': [Ramirez-Ruiz et al. 2005): GRB050525A with SN2005nc (|Della Valle et aLll2006a[) : 
and GRB060218 with SN2006aj (jCampana et al.i ,200 6). In addition, LGRBs in general are associated with low 
meta l licity star-forniing regions in their host galaxies fProchaska et all l2004l : iGorosabel et al.ll2005t [Sollerman et alJ 
2005*: 'Fruchter et al."2006l: iFvnbo et al.ll2006[ ) suggesting a link to massive stars ([Bloom. Kulkarni fc Diorgovskill2002l : 
Wooslcv & Bloom 2006). 

The observational as sociation with SNe ha s given support to the idea that LGRBs are produced in the collapse of 
massive stellar cores. iMacFadven fc WooslevI ([1999f ) developed this in great detail in the coUapsar scenario, assuming 
a black hole (BH) is formed at the center and evaluating the potential to produce a GRB. The infalling stellar material 
that is not ejected in the usual SN forms a massive, dense and hot accretion disk which releases the necessary energy to 
power the burst. They distinguished between Type I and Type II coUapsars: the former entails the direct collapse of 
the iron core to a BH, while in the latter initially a proto-neutron star lies at the center of the star, and later collapses 
to a BH once it has accreted a sufficient amount of mass. A GRB may be produced even if t here is no BH at th e 
center in a magnetically dominated explosion during the proto-neutron star phase of the collapse ([Dessart et al.ll20d8[) . 
Here, we will focus on Type I collapsars. The absence of H lines in the observed GRB/SN spectra implies that the 
progenitor has somehow lost its envelope, perha ps through interaction with a binary companion and/or st rong winds, 
and is thus essentially a Wolf-Rayet (WR) star (Izzard. R amirez-Ruiz fc Toutj|20o4 ICantiello et al.ll2007[) . 

In the standard version of the coUapsar, the GRB/SN link is a quite natural consequence, and observed associations 
have thus provided a strong motivation in this sense to study it further. Typical LGRBs have z ~ 1 — 4, which is 
too far for a SN counterpart to be detected, and indeed, the cases in which there is such a correspondence all lie 
at fa i rly low redshift Two recent events however GRB060505 and GRB 060614 ([Fvnbo et al.l [20061 : iGehrels et all 
l2006t [Gal Yam et al.ll2006[ : [Delia Valle et"aI] l2006bV have pointed perhaps to a different situation in which the GRB 
is not accompanied by a SN. Their low redshift allowed deep searches which convincingly ruled out the presence of 
an underlying type Ic stellar explosion of the kind seen in SN1998bw and SN2003dh, an d particularly, the host of 
GRB060505 is a star- forming galaxy, similar to that of typical LGRBs ([Fvnbo et al.[[200"6[ ). We explore here how the 
initial angular momentum distribution in the star, and possibly the vigor of angular momentum transport within the 
centrifugally supported accretion disk, may furnish a key ingredient to understand this behavior. 

An important point is that despite considerable effort (see .Wooslev fc Heger (2006)), the stellar rotation rate in 
pre-SN cores is not fully determined. It can sensitively depend on evolutionary details, such as mas s loss on the 
main sequence an d the stellar magnetic field, which both lead to important spin down in later stage s ([SpruitI [20021: 
Heger et al.l [20051) . Binary interactions ma y als o affect the rotati o n rat e through tidal interactions ([Detmers et al.l 



20081 ). We note that lYoon fc Langed ([20051 ). and [Wooslev fc Heged ([20061 ). have identified a channel for massive stars 



in which mass and angular momentum losses are greatly reduced by complete mixing in the main sequence and the 
consequent absence of a giant phase. 

What we do know is that the specific angular momentum, J, needed to form a disk around the BH is at least the 
value of the innermost stable circular orbit Risco {Risco — 3rg, where rg = 2GMbyi/c^ is the Schwarzschild radius) 
times the velocity of the particles (which are very close to the speed of light). So, a first estimate for the critical J 
would be Jcrit — RiscoC — 3rgC = 6GM/c ~ 10^^ cm^ s~^ for a one solar mass BH (the true critical value is Jcrit 
2rgC. see Sextion 3.1). Most studies o f collapsars and neutrino-dominated accretion flow s (MacFadvcn & Woosle 
19991: [Popham. Wooslev fc Frved 11999'; 'Heger et al.[ [20001: [Naravan. Piran fc Kumail [2001t [Proga fc Begelman .200g 



Proga. MacFadven. Armitage fc Bcgclman 2003; L ee. Ramirez-Ruiz fc Pagel2005l ) have considered angular momentum 



distributions that are well above this limit and essentially guarantee the formation of a centrifugally supported accretion 
disk (in stellar core collapse calculations as well as disk evolution studies). In practice, there is also something akin to 
a maximum value of angular momentum for the coUapsar model to work, since for very rapi d rotation the accretion 
disk f orms at large radii an d the binding energy cannot be effectively dissipated as neutrinos ([MacFadven fc Wooslev[ 
[19991 : [Lee fc Ramirez- Rui3 [2006). Clearly in the limit of slow rotation the accretion disk will form very near the BH, 
and so general relativity (GR) will play an important role in its evolution. 

Our general objective is to explore the morphology, energy release, and observable signatures of hypercritical ac- 
cretion flows from massive stellar core collapse for a range of angular momentum values covering the transition from 
centrifugally supported disks to near radial inflow (lying below those usually considered in coUapsar calculations) in 
order to better u nderstand the possible production of LGRBs as well as their putative progenitors. We improve upon 
previous studies ([MacFadven fc Wooslev[ [19991 : [Proga. MacFadven. Armitage fc Begelman|[2003[ : [Lee fc Ramirez-RuizI 
[2006) in signiflcant ways by considering more detailed thermodynamics in the equation of sta te, an improved treatmen t 
of neutrinos, and through realistic initial conditions taken from evolutionary calculations of Woosl ev fc Heged ((2006). 

We first describe the input physics and numerical setup in Section [21 followed by our results in Section (31 Prospects 
for GRB production and the link between SN and GRBs are given in Section [D 



2. INPUT PHYSICS 

2.1. Equation of State and cooling processes 

We use a detailed equation of state ([Lee. Ramirez-Ruiz fc Pag3[2005[ ) where the total pressure, P, contains contri- 
butions from an ideal gas of a particles and free nucleons in nuclear statistical equilibrium (NSE), -Pgas, blackbody 
radiation, Pj-ad (the optical depth to photons in the gas is such that they are fully trapped), neutrinos, Pj/, and 



relativistic electron/positron pairs of arbitrary degeneracy, Pf,± ()Blinnikov. Dunina-Barkovskava &: Nadvozhinlll996f) . 
We allow for ncutronization and a correspondingly variable electron fraction Yq by requiring charge neutrality and 
equilibrium in weak interactions, depending on wheth er the fluid is optically thick or thin to its own neutrino emission 
(|Beloborodovll200a iLee. Ramirez Ruiz fc Pag3l2005[ ) . 

The physical conditions in the accretion flow within the collapsing stellar core are such that the temperature and 
density is T ~ 10^ — 10"'^° K, and p ~ 10^ — 10"'^° g cm~^ respectively. Thus the reaction rates that will dominate the 
neutrino emissivities, q^, are capture onto free nucleons {(jcap) and annihilation {qann): for which we use the 
tables of iLanganke fc Martmez-Pinedol (|2P01) and the fltting functions of Utoh et all ()1996f ). respectively ^. Finally, 
photodis integration and synthesis o f a pa rticles can also cool or heat the gas, and is correspondingly accounted for 
(see also lLee. Ramirez- Ruiz fc Pagd (|2005[ )). 

If the fluid is optically thin, the total cooling is simply the sum of the emissivities described above. At a finite 
optical depth, we may split it into scattering and absorption components as t^, = Tgcat + '''abs, where the absorption 
term is due to the inverse reaction of capture onto protons or neutrons. Tabs-cap, and lyV annihilation, Tabs-ann 
pi Matteo. Perna fc Naravan|[2003 ): 

Tabs Tabs— cap ^" Tabs-ann: (-1) 

with 

_ 2 qcap H 

Tabs-cap — „ rriA \^) 

7 CTSB J 

2 qann H 

Tabs-ann — ■= \'J ) 

I OSB J 

H being the pressure scale height in the disk and crge the Stefan— Boltzmann constant. 
The contribution from scattering off free nucleons is 



a 

m„ 

where y„ and Yp are the neutron and proton fractions, Cg^p = (1 -I- 5a^)/24, Cg^n = [4(C^ — 1)^ -I- 5a^]/24, Cy — 
0.5 + 2 sin^ 6I„, sin^ 9^1 « 0.23, a = 1.76 x 10"'''* cm^, and a = 1.25 (jShapiro fc Teukolskvl [19830 and the rest of the 
symbols have their usual meanings. 

To compute the i nternal energy, cooli ng rate, and p ressure due to neutrinos we used a two -stream approximation 
(jPopham fc Naravan,1995. : .Pi Matteo. Perna fc Naravan 2002 : iJaniuk. Yuan. Perna fc Pi MatteQ.2007,) : 

, , / ^ 7^HkBT)'' f (Tabs+Tscat)/2 + l/V3 \ .3 

El, (t) = 3Py (t) = - — V,^ ,„ ^ r=r^ erg cm , (5) 

' ^ ' 8 • 15(fic)3 (rabs + T,eat)/2 + 1/ VS + l/3Tabs / 

7/ 4asBrV3 _3 _i ... 

Qu [t] — -A 1= — erg cm s , (6) 

8 V (Tabs + Tscat)/2 + 1/V3 + l/3Tabsy Vcm; 

and 

1 fTj{^{kBTf\ ( (Tabs+Tscat)/2+l/V3 \ 

^ 3 ^8-15(MV lvKbs+Tscat)/2+l/V3+l/3Tabsy' ■ ^ ' 

The total neutrino luminosity (in erg s~^) is then calculated through 

L(t) = / g,(T) dV erg s-i. (8) 



T,eat = 13.8 (CpFp + C,,„r„) — {k^T/m.c^f pH, (4) 



2.2. Initial setup and numerical method 

Our initial data were taken from the one-dimensional pre-SN calculations of lWooslev fc Hegeil (120061) . Specifically we 
considered model 16TI, a rapidly rotating {v^ot = 390 km s^^ at the equator), 16 Mq WR star of low metallicity (1% 
solar), low mass loss (2 Mq loss at the end of the evolution), with an iron core surrounded by silicon and oxygen, neon 
and carbon shells. The density, temperature, and radial velocity distributions as functions of the spherical radius, 
i?, were mapped onto a two-dimensional configuration assuming spherical symmetry. We further suppose that the 
low angular momentum iron core (with mass Mpe = 1.6Mq) will promptly collapse to a BH (possibly producing a 
Type I collapsar, see Section [T]) and thus condense all of this matter onto the origin in a point mass at the start of the 
simulation, when the surrounding envelope begins its infall. 

The calculatio ns are performed with a two-dimensional Smooth Partic le Hydrodynamics (SPH) code 
(IMonaghanI |1992[) in cvlindrica l coordinates (r, z) with azimuthal symmetry (|Lee fc Ramirez-RuizI l2002l 120061 : 
iLee, Ramirez-Ruiz fc Pagil2005l ) and a reflecting boundary along the polar axis (no reflection about the 2 = plane 



^ Neutrino emission through plasmon decays and nucleon— nucleon bremsstrahlung was also computed and found to be insignificant 
compared to capture and annihilation. 
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Fig. 1. — Density, temperature and radial velocity profiles in the pre-SN star used as an initial condition (model 16TI of lWooslev &: Hegeil 
l|2006l )'). The solid vertical line in each panel marks the outer computational boundary used in the simulations. 



is assumed). The inner and outer boundaries were set at spherical radii Rin = 2 x 10® cm and i?out — 2 x 10^ cm, 
respectively, with free outflow conditions imposed during the simulation. Any material crossing the inner boundary is 
assumed to be accreted by the BH, whose mass is correspondingly updated. Matter lost through the outer boundary 
is not followed. We solve d the continuity, moinentum , and energy equations, includ ing the terms arising froin the full 
viscous stress tensor tr^ (|Lee fc Ramirez-Ruizjl20d^ . and took the a prescription (IShakura fc Sunvaevlll973h for the 
coefficient of viscosity, Tyv = apcs H, where Cg is the local sound speed and H = Cs/^k is the pressure scale height 
(Qk is the local Keplerian orbital angular frequency). Doing the calculation in two dimensions allows for good spatial 
resolution, a considerable simulated interval and a solid discussion of angular momentum and viscosity effects. The 
initial profile is reproduced with a Monte Carlo accept/reject procedure, typically with 1 — 2 x 10^ SPH particles in 
the computational domain. The spatial resolution is intrinsically adaptive in this numerical scheme and varies greatly, 
increasing at small radii (or equivalently, high densities) with a typical smoothing length of few x 100 m in the inner 
disks formed. 

The gravitational potential of the BH is computed with the pseudo-Newtonian expression of iPaczvhski fc Wiital 
(|T980l) : 

<^ = -^, (9) 



which approximates GR effects in the inner regions for a nonrotating (Schwarzschild) BH. In particular, it reproduces 
the position of the innermost stable circular orbit at Visco = 3rg, as well as the marginally bound orbit at = 2rg. 

The mass in the collapsing envelope is comparable to that of the BH, so one should consider self-gravity. We assume 
that the mass distribution remains roughly spherical, so a mass element at spherical radius i?o is affected only by the 
matter distribution at radii R < Rq as if it were concentrated at the center of the star. This, while strictly valid only 
for configurations with spherical symmetry, is nevertheless a good approximation for the present set of calculations, 
since deviations in the mass distribution indeed remain small. Figure [T] shows the initial density, temperature, and 



radial velocity profiles in the pre-SN star as functions of spherical radius, taken from IWooslev fc Hegerl ()2006[ ). 



3. RESULTS 
3.1. Expectations 

Consider a nonrotating BH surrounded by rotating material distributed in spherical symmetry with vanishing pres- 
sure. For a distribution of specific angular momentum increasing monotonically with respect to the polar angle, 
flowlines in the polar regions will be qualitatively different than those near the equator, and can be divided into three 
types: (1) those with very little angular momentum, moving nearly radially into the BH; (2) those with high angular 
momentum which can find a point where centrifugal support balances the gravitational field (the circularization ra- 
dius); and (3) those that would accrete onto the central mass but cross the equatorial plane, z = 0, before doing so. 
Recall that in GR capture orbits exist even for nonvanishing angular momentum, and the critical value defining this 
class of solutions is Jcrit — 2.0rgC. Here we will consider J > Jcrit as high ang ular momentum, and J ^ JctH. a s low. 
This general situation has been cons idered before in the context of LMXBs by iBeloborodov fc lUarionovl (|2001[ ) , and 
for accretion modes in coUapsars by 'Lee & Ramirez-Ruiz ("2006") (an analytical derivation of low angular momentum 
flow lines in GR has been carried out by Mcndoza ct al. (2009)). 

Encountering the centrifugal barrier or a flowline from the opposite hemisphere can lead to a shock (and thus 
hydrodynamical effects) in which the kinetic energy is efficiently converted into thermal energy. If the material does 
not have enough angular momentum to remain in orbit, it will fall onto the BH due to GR effects, even in the absence 
of any angular momentum transport mechanism through a fast, nearly inviscid disk. For high angular momentum, a 
disk (with scale height H and typical radius r) forms in the equatorial plane, and if cooling occurs some of this energy 
will be lost from the system through photons or neutrinos, depending on the physical conditions. For efficient cooling 
the disk will be nearly isothermal and geometrically thin, with _ff <C In the opposite adiabatic regime it can be quite 
thick, and H ^ r. Energy dissipation will in general circularize motions at a radius given by the local value of angular 
momentum, and if any of this is subsequently removed (through any mechanism), or it is already below the critical 
value, the material inside the disk may move radially and possibly be accreted by the BH. The size and geometry of 
the newly formed disk thus reflects: (1) the angular momentum distribution, fixing the point where the centrifugal 
barrier is encountered; (2) the mass accretion rate, giving the total gravitational energy which can be converted into 
thermal energy; (3) the cooling rate, determining the disk geometry; (4) the value of the a parameter, responsible for 
angular momentum transport. 

In gene ral, previous collapsar studies or neutrino-dominated accretion flow studies have considered cases where 
J ^ Jcrit f MacFadven fc Wooslcv 1999: Poph am. Wooslev fc Fr ver 1999: Heger ct al. 2000; N aravan. Piran fc Kumail 
1200 It iProga fc Bcgelma n 2003; Proga. Mac Fadven. Armitaee fc Begelman 2003). Only recently has the low angular 
momentum regime in the collapsar context been addressed ( Lee fc Ramirez- Ruid 120061 ) . showing the presence of the 
small dwarf disk as described above. For J > 2.0rgC, a thick disk forms, supplying th e bulk of the neutrino luminosity 
from a shocked toroidal region around the BH. More recentlv. I Janiuk fc Progal (j2008l ) considered the minimum angular 
momentum required for disk formation, taking into account that the accretion of mass by the central object raises the 
angular momentum threshold as the stellar core collapses. 

We explored the dynamical evolution for different values of the a viscosity parameter (see Table 1), and more 
importantly, different values of angular momentum, which were always separated into radial, R, and polar angle, 9 
(measured form the rotation axis), variations as J = J{R,d) ~ Jb.{R) Je{0)- For the radial component, J^{R), we 
considered constant values with radius (see Table 1). We will consider more realistic distributions of Jn{R) in future 
work (D. Lopez-Camara et al. 2008, in preparation). In all cases the initial distribution for J (9) corresponded to rigid 
body rotation on shells, J{9) = sin^ 9. 
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3.2. Global properties and flow morphology 

For J(i?) = 2.0rgC, the flow is essentially at the critical value, and remains nearly radial. Even though some 
compression occurs in the equatorial region the centrifugal barrier is absent and no shocks are formed. Moreover, 
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Fig. 2. — Velocity field and density map (in g cm ^) for ( Jq, a) = (2rgC, 0.1) at t = 0.2 s. The largest vector corresponds to 8 X lO'^ 
cm . 

since J(i?) is constant, the initial velocity field and density map (shown in Figure [2]) do not vary substantially as the 
simulation progresses. For an inflow in strict free fall, conservation of mass and energy give a density profile 

-3/2 / Mw ^ "^^^ 



P = 



M 



1.2 X 10^ 



M 

0.5 Mm s- 



-1 ) (lO^cm) 



^BH 

1.7 Mp 



g cm 



(10) 



where M is the mass accretion rate and Mbh the BH mass. The solution is plotted as a dotted line in Figure |3] with 
the above scalings set to unity, clearly showing the excellent agreement between the full simulation and this simple 
estimate along the equator. The energy release, however, is not negligible, as will be seen below in Section [331 As 
long as the equatorial angular momentum is below this critical value, the solution is insensitive to the actual value of 
the viscosity parameter, a, that is chosen (we computed solutions for a in the range — 10~^). 

For values of J(i?) above the critical value the outcome changes drastically. The scale of the disk correlates with Jq 
so we will limit the discussion to the case with J(i?) = S.OrgC and a — 0.10. An initial disk forms close to the equator 
after 0.1 s, the free-fall time from the boundary of the Fe core, when the innermost regions of the envelope approach 
their circularization radius Vc ~ J§ /GMbh- The transformation of kinetic into thermal energy at the centrifugal 
barrier produces a hot torus around the BH. Figure H] shows the velocity field and a density map at t = 0.2 s for 
(Jo, a) = (3rgC, 0.1). Note how in the velocity field the fiow lines in the polar regions are still largely radial, and 
an asymmetry is apparent, with the polar shock lying substantially closer to the black hole than in the equator. A 
large-scale meridional circulation with two prominent eddies is also apparent. 

The temperature and density rise rapidly behind the shock front, releasing a substantial amount of energy and 
photodisintegrating He. It is here where neutrino emission, finite optical depth effects and neutronization can be- 
come important. A substantial fraction of the released energy will likely be directed to the polar regions because of 
geometrical effects, and may give rise to a fireball capable of producing a GRB (see Section [375]) . 

Outside the shock front, the solution is essentially the same as for low angular momentum (see Figure [3]), with a 
jump at r « Tc. For angular momentum Jq = Ar^c, with A > 2 being a constant, the circularization radius for 
equatorial matter is Vc w '^•A^rg and the free-fall velocity at Tc is Vff = c/2^/^A. If the kinetic energy of infall is entirely 
transformed into thermal energy, one would naively estimate the latter as 

2 

kT^^^10iA/3)-'MeY. (11) 

This is clearly an overestimate of the temperature, since the velocity is not entirely radial but has a substantial 
azimuthal component, v^, by the time it reaches the circularization radius, rc- It does nevertheless give a useful 




Fig. 3. — Equatorial {z = 0) density, temperature, electron fraction and nucleon mass fraction distribution for {Jo, a) = (3rgC, 0.1) (solid 
line) and ( Jq, a) = (2rgC, 0.1) (dashed line) at t = 0.2 s. The power law given by a thin solid line in the density panel is the solution given 
in Equation (10). 

estimate (and upper bound) on the postshock temperature, as can be seen from Figure [3l As we shall see, it is also 
high enough to guarantee the complete photodisintegration of the infalling nuclei into their constituent neutrons and 
protons. 

In the initially small disk that forms at the centrifugal barrier, if a 7^ transport processes transform the constant 
distribution of angular momentum into a nearly Keplerian one where differential rotation is important. For the 
Newtonian case the orbital frequency is cx r^^^^ and at a few gravitational radii this is not a bad approximation 
even in the pseudo-GR PW potential. In a way this partially erases the memory of the initial distribution: as long 
as the rotation rate is above a certain threshold, the inner disk will quickly converge to a ccntrifugally supported 
structure with strong differential rotation, independently of the details of the initial angular momentum profile. The 
associated shear produces dissipation and pumps mechanical energy into thermal energy. The rotational structure of 
the inner disk depends on whether a is finite or not. If a 7^ the actual value of the specific angular momentum 
can grow beyond the initial equatorial value Jq through transport. For the inviscid regime, Q oc r~^/^ (and hence 
J cx r^/^) only as long as this implies J < Jo- At greater radii the equatorial flow maintains constant specific angular 
momentum. 

Approximately 0.2 seconds after the formation of the initial disk, a shock begins to propagate outward, initially 
moving at u w 2 x 10^ cm s~^ in the equatorial regions. The postshock gas pushes out against the infalling envelope 
mainly because of viscous heating for a = 0.1, and is aided by a combination of neutrino heating and He synthesis ^. 
The calculations were stopped at i ~ 0.5 s, when the shock reached the outer boundary. Not only is the numerical 
solution no longer self-consistent, but by then our resolution had decreased notably, since accretion onto the BH entails 
a loss of particles, and thus resolution. We will address this issue in an adaptive form in future work. For low viscosity 
the inner disk is denser due to the reduced transport efficiency and the outward motion of the shock is only slightly 
delayed. In the inviscid limit there is no associated heating but the piling up of material due to the lack of transport 
also induces outward motion of the shock front after a comparable delay of ~ 0.2 s. 

In general, the inner disk can remain thick despite the presence of cooling because of the continous infall of material 
and, in some cases, because the optical depth to neutrinos is large enough to keep the internal energy from immediately 



^ In a test calculation with the viscous heating terms switched off, the shock had not started to move outward after 0.5 seconds. 
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Fig. 4. — Velocity field and density map (in g cm ^) for (Jo, a) = (3rgC, 0.1) at t = 0.2 s. The largest vector corresponds to v 8 X 10^ 
cm for densities lower than 5 X 10* g cm""^. At higher values the velocity vectors are reduced by a factor of two. The location of the 
accretion shock is clearly seen at r{z = 0) ~ 2 X 10^ cm. 

escaping. As long as the angular momentum of the infalling gas is greater than the critical value, the morphology is 
similar, with a greater radial extent for higher Jq. 

The breaking of spherical symmetry is in principle a problem for the computation of self-gravity, described in 
Section [2?2l We thus checked how far the actual flow departs from spherical symmetry in terms of the mass distribution 
as a function of the polar angle 6, and find that the scatter in the radially integrated mass along cones with constant 
6 is small enough (10%) to be ignored in a first treatment. 

3.3. The importance of cooling 

To highlight and better understand the importance that the proper computation of cooling has on the global 
properties of the solution, we have calculated the evolution of the flow in two simplified cases. Since neutrinos are the 
only means other than advection onto the BH through which the cas can cool, and thus move lower in the gravitational 
potential well, there are two limits in this respect: adiabatic inflow, in which no cooling occurs, and isothermal flow, in 
which on the contrary, it is extremely efficient. The true solution must lie somewhere between these two extremes, and 
it is instructive to know which it resembles the most. We computed these, in one case by eliminating the cooling terms 
in the energy equation and thus impeding the outward flow of energy through neutrinos, and in the other by using and 
ideal gas equation of state with P = (7 — where 7 = 1.01 and the same initial conditions, thus mimicking the 

isothermal case where 7 = 1 and compressibility is very high. The velocity fields and the corresponding density maps 
(also for Jo = 3rgC and a — 0.1 are shown in Figures O and El Equatorial density profiles for the various cases are shown 
in Figure [71 and can be compared with those in Figure [2]). The isothermal fiow looks qualitatively similar to the low 
angular momentum case, but the density is 2 orders of magnitude larger. There is little temporal evolution once the 
centrifugal barrier is reached, and the solution is quasi-stationary. On the other hand, in the adiabatic solution as soon 
as the infalling gas reaches the centrifugal barrier, the shock bounces rapidly outward and the subsequent expansion 
produces a strong flow reversal, sweeping the inner envelope outward and beyond the outer boundary (compare with 
Figure [Hand note that the spatial scale is 5 times larger). At the instant shown in Figure [7| the density has decreased 
considerably, even beyond the solution where neutrino emission was included. The fact that the accretion disk remains 
geometrically thick in the calculation with neutrino cooling indicates that it is not extremely efficient, resembling the 
adiabatic solution qualitatively, but significant enough to avoid such a prompt outfiow. Given that the emissivities 
are sensitive functions of temperature and composition, an accurate expression for these and at least an approximate 
treatment of neutrino optical depths is clearly a crucial ingredient in the evolution. Note that since the optical depth 
to neutrinos is never extremely large, once emitted they are relatively free to escape and are not advected with the 
flow. 
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Fig. 5. — Velocity field and density map (in g cm~'^) for isothermal flow with Jq = 3rgC and a = 0.1 at t = 0.2 s. 

3.4. Thermodynamics inside the disk 

The origin of the energy release and the associated neutrino luminosity can be easily understood by examining the 
various emission mechanisms separately. Figure |H] shows a snapshot at t = 0.2 s of the neutrino emissivities for pair 
annihilation {(jann) and capture onto free nucleons {qcap) as a function of the cylindrical radius r in the equatorial 
plane (z=0), for {Jo, a) = (3rgC,0.1). 

The highest energy release occurs very close to the BH: at r < 2 x 10^ cm both qcap and qann rise at least 3 orders of 
magnitude and release up to 10^^ erg cm~'^. The figure also shows the corresponding optical depth contributions 
for both processes as well as coherent scattering (rgcat) for (Jo, a) = (3rgC, 0.1) a,t t = 0.2 s. The inner regions 
(r < 2 X 10^ cm) can become somewhat opaque while at large radii the fluid remains optically thin. Electron/positron 
annihilation dominate the optical depth in the dense inner region over captures onto free nucleons and coherent 
scattering. Neutrinos are not entirely free to escape and corrections due to finite optical depth both in the pressure 
and luminosity (Eq. [6]) must be taken into account. A clue to the disk's structure can be inferred from the comparison 
of the local coohng (or Kelvin-Helmholtz) time, tcooi ~ E/q^ to the dynamical time, tdyn ~ {GMBii/r^)~^^'^ , plotted 
in Figure [HI Despite the high accretion rates, even in the inner regions the cooling time is much longer than the orbital 
period ^. This is in agreement with the previous comparison to the isothermal and adiabatic limits in Section 13.31 
where we already found that qualitatively at least, the overall picture is indicative of inefficient cooling. 

Regarding the thermodynamical properties of the flow, we will focus on ( Jg, a) = (3rgC, 0.1) and ( Jq, a) = (2rgC, 0.1) 
(the remaining cases are quite similar). Equatorial (z — 0) profiles for density, temperature, electron fraction and 
nucleon mass fraction are shown in Figure [3] for both runs. The outwardly moving accretion shock front in the high 
angular momentum case is visible at r « 2 x 10^ cm. Once in the postshock region, compression and the rise in 
temperature fully break up a particles into their constituent neutrons and protons. The consequent cooling through 
captures onto neutrons and protons allows the gas to reach densities high enough that neutronization takes place, 
lowering to minimum values close to 0.2 (the higher densities also make the fluid become more degenerate). 

For the calculation with J{R) = Jmt no centrifugally supported disk forms, although compression is important in 
the flow just before the gas falls through the inner boundary. The equatorial profiles (shown in the previous figures 
along with those for Jq = S.OrgC) show little evolution, and this case is basically in the inviscid limit, where varying 
a makes little or no difference on the outcome. The gas is so close to being in free fall that the much longer viscous 
timescale becomes irrelevant, and runs with different values of a yielded the same solution. 

Contrary to the cases with J(-R) > Jcrit, where the disk was in some cases opaque, the maximum optical depth is 




At large radii it is appropriately even larger, since the stellar envelope is essentially in hydrostatic equilibrium. 
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Fig. 6. — Velocity field and density map (in g cm for adiabatic flow with Jo = 3rgC and a = 0.1 at t = 0.2 s. Note the different spatial 
scale compared to the previous figures. The largest velocity vectors correspond to lO^cm for densities lower than 5 X lO'^g cm~'^. At 
higher values the vectors are reduced by a factor of 2. 
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Fig. 7. — Equatorial {z = 0) density distributions at t = 0.2 s for the simulations carried out in the adiabatic (dotted) and isothermal 
(dashed) limits, both for (Jo,a) = (3rgC, 0.1). The thick solid line is the simulation with full microphysical treatment and (Jo, a) = 
(SrgcO.l). 



now at most r = 0.01, dominated also by pair annihilation (the temperature is kT ~ 4 MeV) and the emissivities 
are substantially lower. Compression and the associated rise in temperature are strong enough to dissociate He 
almost completely in the inner regions, but no significant neutronization occurs because the maximum density is only 
~ 10^ g cm~^. Very close to the BH, e"^ captures produce slight neutronization and the electron fraction drops down 
to Ye = 0.45. The pairs are at the threshold of degeneracy, with /ic/fceT « 1.5, where fie is the chemical potential 
of the electrons. 




Fig. 8. — Upper panel: neutrino emissivities qann (solid line) and qcap (dashed line) for {Jo, a) = (3rgC, 0.1) at t = 0.2 s. Lower panel: 
the corresponding optical depth contributions from Tabs_ann (solid line), Tabs— cap (dashed line), Tscat (dotted line) for {Jo,a) = (3rgC, 0.1) 
at t = 0.2 s. 



The possibility of this transition at high accretion rat es (of o rder a few tenths of a solar mass per second here) 
was hinted at in the simpler calculations of Lee & Ramir ez-Ruid ([2006), but the equation of state used there did not 
permit neutronization or degeneracy effects to play a role ^. It is clear from the present set of calculations that such 
approximations are valid only for quite lower accretion rates (a few hundredths of a solar mass per second). This 
limitation is valid also for other collapsar studies, where neither finite degeneracy or neutronization have been fully 
considered. Two important potential consequences in this context stand out. 

First, the energy of a given neutrino depends on the physical process responsible for its creation. Those arising 
from pair annihilation (a thermal process) have characteristic energies of the order E^^nn ~ 4fcr, while those due to 
capture onto free nucleons (a weak interaction) have energies of the order ii^cap ~ E-p, where Ep ~ 9(pioi^)^^'^ MeV is 
the Fermi energy of electrons/positrons (pio — p/10^° g cm~^). Table 2 shows the values for pio, Tio, Ye, Eann, and 
Ecap for our calculations (in the inner regions where the energy release is largest, i.e., where p — Pmax, T — Tmax, and 
Y^ = Y^,min)- As -Eann and i?cap Vary with respect to Jg and a, the transition to a different thermodynamic regime will 
thus modify the emergent neutrino spectrum. This in turn affects the global energy release, since at a given neutrino 
luminosity, ij^, the efficiency for neutrino annihilation scales with < E^ >, and thus more energetic neutrinos will be 
more efficiently converted into a relativistic pair plasma. 

Second, the morphology of the flow may be influenced by the available energy sinks (e.g., advection, neutrino losses) 
allowed in a calculation. Dissipation in the disk through viscosity adds to the internal energy reservoir. If it is allowed 
to escape with moderate efficiency the vertical scale height will be limited. Otherwise the d isk will tend to expand 
in res ponse, or possibly drive winds from its surface. In the original collapsar calculation of iMacFadven fc WooslevI 
()1999r ). strong winds driven from the surface of the disk at small radii were reported for one calculation, with a = 0.10. 
Calculations with a smaller dissipation rate did not exhibit this feature. For now, our computations do not span such 



^ They assumed that hot pairs were abundant, producing a pressure contribution oc T^, and that Ti, <C 1. 
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Figure |4] 
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Neutrino Energetics 
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a long time interval, but the essential morphological features (expanding outer accretion shock, hot torus, and dense 
inner disk) are well established by 0.4 s. We intend to further explore these issues in a second set of simulations 
with varying angular momentum distributions. Any potential outflows emanating from the inner regions could have a 
significant impact on the observable signature of the stellar collapse, and so this deserves careful consideration. 

3.5. Accretion rate, neutrino luminosity, and energy conversion efficiency. 

The mass accretion rate Mbh (t) was computed at the inner boundary R — Rin = 20 km and is shown in Figure [TO] 
for three runs with the same value of a = 0.10, and different angular momentum Jq. The corresponding neutrino 
luminosity is plotted in the top panel of Figure [TTl The initial delay of about 0.1 seconds represents the infall time 
from the outer boundary of the iron core. Thereafter the accretion rate and the luminosity rise rapidly, reaching 
approximately 0.6 solar masses per second and 10^^ — 10^^ erg respectively. The inner disk is responsible for 
most of the energy release, and its formation and steady configuration for the next 0.1-0.2 seconds produce a nearly 
constant accretion rate and luminosity. The efficiency L^/Mbhc^ during this period is approximately 0.01 and 0.001, 
for high and low angular momentum, respectively, reflecting the increased importance of advection when the flow is 
quasi-radial. Note that despite the absence of a centrifugal barrier and a dense disk, the configuration with low angular 
momentum is capable of significant energy release {L^ ~ 10^^ erg s~^). 

The launching of the accretion shock after a few tenths of a second leads to a drop in both the accretion rate and 
the luminosity, while maintaining a nearly constant efficiency. For Jq/ (rgc) = 2.5,3.0 the outwardly propagating 
shock eventually perturbs the mass fiux to the inner disk before reaching the outer boundary used in the calculation. 
The higher mass accretion rate at low angular momentum (below the critical value) is due to GR effects but does 
not initially translate into a higher luminosity because the material falls nearly radially into the BH and is unable to 





Fig. 10. — Mass accretion rate onto the BH (Mbh) for a fixed a with different angular momentum values: (Jo, a) = (2rgC, 0.1) (red solid 
line), (Jo,o) = (2.5rgC, 0.1) (green dashed line), and (Jo,o) = (3rgC, 0.1) (blue dotted line). 

radiate efficiently. 

Neutrino luminosities for a fixed value of angular momentum ( Jq = S.Orgc) and varying efficiency of viscous transport 
are shown in the bottom panel of Figure [TT] A moderate to null viscosity is in this early stage the most efficient at 
producing neutrinos in spite of the lower dissipation rate since very efficient angular momentum transport drains the 
inner disk much too rapidl y. 

Previous studies (jPi Matteo. Perna fc NaravanI |2002| : iLee. Ramirez- Ruiz fc Pagil2004 120051 ) have shown that at 
very high accretion rates, greater than approximately one solar mass per second, the neutrino luminosity saturates 
and t he accretion efficiency decreases, inhibiting the driving of winds from the disk (|Di Matteo. Perna fc NaravanI 
I2OO2D . For the cases computed here, regardless of the actual value of angular momentum considered, the flow never 
becomes strongly opaque to its own neutrino emission in the sense that the emitted energy is able to escape, thus 
avoiding this limitation. Resolving this issue in more detail requires exploring a wider range of progenitors and longer 
timescales than those considered here and is beyond the scope of the present paper. 



4. DISCUSSION AND CONCLUSIONS 
4.1. Limitations and comparison to other work 

As with all numerical work, the choices made in carrying out the simulations reflect intentions and biases, and the 
current investigation lacks in several aspects. For example, being a two-dimensional calculation, self-gravity has been 
considered only in an approximate manner, and stability issues that relate to this or are intrinsically three-dimensional 
(such as spiral arms) are ignored. 

Magnetic flelds have not been included and are in all likelihood important in several aspects, two of which deserve 
special mention. The first is simply that magnetic fields are thought to be at the origin of the magneto rotational 
instability (MRI) (Balbus & Hawlcy 1998) responsible for angular momentum transport and dissipation, which we 
have considered through the a parameter. The MRI likely operates in coUapsars, although the details of how it does 
so are likely to be different than what we can infer about it from theory and observation of accretion disks in CVs and 
X-ray binaries in a more leisurely state of affairs. The second point is relate d to the importance of magnetic pr essure 
in possible outflows emanating from a disk. Detailed GR MHD calculations (IPe Vill iers. Hawlev fc Krolikll2005f l show 
that a flow starting from an equilibrium torus develops the MRI and produces an inner dense disk, a funnel wall, and 
a corona around the central BH. Gas pressure dominates largely in the first, so our hydrodynamical solution to the 
inner disk is likely to be a realistic approximation, while magnetic pressure is relatively more important in the other 
two, a nd being polar structures th ey are important for the driving of outflows and energy release in the context of 
GRBs (jMcKinnev fc Naravan|[2b07[ ) . For very rapidly rotating BHs. iKrolik. Hawlev fc Hirosd (|2005f ) have shown that 
the magnetic field can actually have a dynamical effect on the mass accretion rate, suppressing it in the inner regions 
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Fig. 11. — Upper panel: Neutrino luminosity for a fixed value of a = 0.10 and different angular momentum values: (Jo, a) = (2rgC, 0.1) 
(red solid line), (Jo,o) = (2.5rgC, 0.1) (green dashed line), and (Jo,o) = (3rgC, 0.1) (blue dotted line). Lower panel: Neutrino luminosity 
Lj/ for a fixed angular momentum Jq = S.OrsC and varying a = 0.00 (red solid line), a = 0.01 (green dashed line), and a = 0.10 (blue 
dotted line). 

(in their case the Kerr parameter was a = 0.998), although a full evaluation of this deserves further study. 

GR, approximated in this work simply by use of the Paczyhski-Wiita potential, is likely to play a role as well. The 
location of the horizon and last stable orbit are functions of the rotation of the BH, lying at rn = 2GM/c^,GM/c^ 
and nsco = QGM/c^jSGM/c^ for a = 0, 1 respectively, and this is important since most of the energy is released 
at small radii. For our purposes, however, they are not likely to be crucial since the relevant radii shift appreciably 
only for very rapidly rotating holes (with a above 0.7 or so), and the progenitor stars do not easily reach such high 
values (model 16TI has a = 0.44). A secondary aspect of GR is related to the efficiency for vT' annihilation and the 
corresponding energy release. Some energy is directly lost to the BH because of the strong gravitational field, while 
focusing increases the annihilation efficiency. Recentlv lBirkl. Alov fc Jankal (|2007T ) have carefully computed the energy 
deposition rates and efficiencies for annihilation in the Kerr metric for various BH spin rates and flow geometries. 
They find that the power output is affected by approximately up to a factor of 2, depending on whether the disk is 
geometrically thin or thick (or even for spherical configurations of the neutrinosphere) . Thus the general energy scale 
is affected, but not crucially so. 



On the other hand, our effort has been directed here at improving the thermodynamical treatment and that of 
neutrinos, as weU as in considering reahstic initial conditions derived from steUar evolution calculations. The choice 
of inner boundary means that the dense, inner disk producing most of the luminosity is well resolved and treated 
appropriately. A final point that deserves improvement is clearly the distribution of angular momentum, since having 
a constant value at all radii is not realistic. We have chosen it for now to gauge the effect of other changes when 
compared to previous work, and provide a guide in further investigation, which will fully explore situations in which 
the radial part is a characteristic function of radius (D. Lopez-Camara et al. 2008, in preparation). 

With this in mind, we may consider the implications for GRB production and GRB/SN associations following core 
collapse and prompt BH formation from this set of calculations. 

4.2. Global energetics 

The fact that we have placed the inner numerical boundary fairly close to the BH (at rjn — 20 km) allows us to 
directly and more realistically c ompute the energy rel e ase (o bviously at the cost of a shorter simulation in time). This 
is similar to what was done by iLee fc Ramirez-Ruid (|2006l ) and so we can directly compare the two sets of results. 
The two most important differences between the two studies are in the use of a better equation of state and neutrino 
treatment here, and in the use of initial models taken from stellar evolution calculations. The neutrino luminosities 
are higher in the present study, partly because of the more realistic physics employed, but also because of the initial 
conditions, which directly result in higher accretion rates and densities. Since the cooling rate from captures scales 
with the density, this fact alone will raise the luminosity. The rates are quite sensitive to the temperature, but the 
rise in the postshock material is mostly due to the infall kinetic energy per unit mass and is fixed by the potential well 
of the central BH. It is clear as well, from the results presented for simplified cases in the adiabatic and isothermal 
limits, and from the computation of the Kelvin-Helmholtz timescale, that correctly accounting for the energy sinks 
is a crucial ingredient if one wishes to estimate the global energy release and flow properties. Finally, the emergent 
neutrino spectrum is a function of the angular momentum of the infalling gas, through the selection of the dominant 
cooling mechanism. This is admittedly a difficult issue to resolve observationally, but may impact upon the energy 
deposition rate. 

The two most important variables that determine the global morphology are the rotation rate, quantified here through 
the distribution of specific angular momentum, and the strength of angular momentum transport, parameterized 
through the prescription of Shakura & Sunyaev. It is interesting to note that for high angular momentum (larger than 
the threshold for disk formation), viscous dissipation is responsible for the production of large-scale flows, namely by 
producing an outward moving shock after a certain delay (a few tenths of a second) . The effect is to perturb the fiow 
into the inner disk and decrease the luminosity. The impact on the flow at large radii remains to be explored in greater 
datail due to the limitation of our outer boundary. 

4.3. Implications for GRB production 

The compression for low angular momentum cases releases a large amount of energy, not more than 1 order of 
magnitude smaller than for the case with a disk, and we believe this to be an important issue: slowly rotating models 
are in principle capable of releasing an amount of energy that could produce GRBs, although admittedly perhaps only 
at the faint end of the luminosity distribution. For a 1% efficiency of conversion into pairs through vv annihilation 
at L,y = 10^^ erg s^^, the models shown here can produce annihilation luminosities of about 10^^ — 10^^ erg s~^, 
depending on the value of Jq. Performing realistic runs for longer times is a priority, but clearly over long timescales, 
having slow rotation is not necessarily a handicap r egarding the energy re l ease. The simple reason is that, as known for 
a long time, and first quantified in this scenario by lLee fc Ramirez-Ruid (|2006l ). it is most efficient in terms of energy 
extraction to have the material release its energy as close as possible to the BH, or equivalentely, as deep as possible 
in the potential well, but not too close so that it is still shocked by the presence of a centrifugal barrier. Magnetic 
mechanisms may also power fully or partially a potential burst, and we note here only that the internal energies are 
high enough, as in the most common hypercritical accretion scenarios, to do so if a fraction of this is transferred to a 
magnetic field (at a level of 10% of equipartition) . 

The obvious implication in terms of the type of progenitor one can consider for GRBs is that single stars at the 
threshold for disk formation through centrifugal support are capable of giving a GRB. Having a substantial amount 
of rotation will obviously guarantee an accretion disk, but perhaps not all cases require such special conditions (as 
for example, torquing of the pre-SN star by a binary companion). The isotropic equivalent energy of GRB060505 is 
log[£'^^iso] — 49.5, clearly within the range obtainable with the models presented here even at low rotation rates. 

4.4. Nucleosynthesis in the outflow and GRB/SN association 

Centrifugally supported coUapsar disks are expected to produce strong winds (|MacFadven fc Wooslevl[l999l ). driven 
quite generically by a combination of viscous, neutrino, and magnetic effects (the first two of which we have explicitly 
included here). iPruet etldl (j2004l ) i n particular, computed the expe cted nucleosynthesis of ^®Ni in these outfiows, 
using the steady state disk models of lPopham . Woosle v fc Frveil ()1999f ) as an initial condition and considering various 
mass accretion and dissipation rates. The essential feature of these models is the computation of the change in entropy 
per baryon, Sb/k in the outfiow, starting from the midplane of the accretion disk and reaching a large distance at 
which the velocity approaches an asymptotic value i^w/c ~ 0.1 — 0.2. From our hydrodynamical determination of 
the disk structure using an improved equation of state and thermodynamics, we computed the equivalent entropies 



in the disk and applied the formalism of iPruet et al.l (|2004f ) to compute the total change in Sh/k (assuming the wind 



Ni, 



reaches similar velocities, unresolved in our simulations spanning 0.5 s). The mass fraction converted to ^^Ni, X 
then depends sensitively on Sb/fc and the ratio (3 — Mq/w^, where Mo is the mass outflow rate and Ww is the terminal 
wind velocity. Once the outflow is launched we can obtain an estimate of directly from the momentum field in 
the simulation, and v^/c is taken as ~ 0.2. We find that for Jq > Jcrit, and independently of the strength of angular 
momentum transport, substantial Ni synthesis occurs (X(^^Ni) « 0.5) in an outflow with M w 0.3 Mq s"^ This 
would imply a total ejection of ~ Mq in ^^Ni in a GRB lasting 10 seconds, in the right range to account for the 
outflows seen, e.g., in SN2003dh. The relatively narrow range in Jq spanned in our calculations does not allow for a 
good extrapolation to higher rotation rates, but clearly the amount of mass involved is significan t, and would likely 
grow substantially at higher Jq where more masive and more radially extended disks are expected (iPruet et al.l (l2004l ) 
estimate several Mq for the massive, centrifugally supported disks taken from lPopham. Wooslev fc Frveil ( 19991 )). A 
configuration in which no significant outfiows occur because the rotation rate is too low (J < Jcrit) would not be 
expected in this scenario to produce any significant ejection of radioactive elements capable of producing a SN-like 
signature. It could, as shown here, still release enough energy through neutrinos or magnetic mechanisms to power a 
classical GRB. The result, if it were to occur at low redshift, could possibly be an event resembling GRB060505. 
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